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Summary. In this contribution we review the theory of integrability of quantum 
systems in one spatial dimension. We introduce the basic concepts such as the 
Yang-Baxter equation, commuting currents, and the algebraic Bethe ansatz. Quite 
extensively we present the treatment of integrable quantum systems at finite tem- 
perature on the basis of a lattice path integral formulation and a suitable transfer 
matrix approach (quantum transfer matrix) . The general method is carried out for 
the seminal model of the spin-1/2 XX Z chain for which thermodynamic properties 
like specific heat, magnetic susceptibility and the finite temperature Drude weight 
of the thermal conductivity are derived. 



1 Introduction 

Integrable quantum chains have continuously attracted attention, because of 
the possibility of obtaining exact data for the spectrum and other physical 
properties despite the truely interacting nature of the spins resp. particles. 
Important examples of these systems are the Heisenberg model [1,2], tJ- [3-5] 
and Hubbard models [6, 7] The computational basis for the work on integrable 
quantum chains is the Bethe ansatz yielding a set of coupled non-linear equa- 
tions for 1-particle wave- numbers (Bethe ansatz roots) . Many studies of the 
Bethe ansatz equations were directed at the ground-states of the considered 
systems and have revealed interesting non-Fermi liquid properties such as 
algebraically decaying correlation functions with non-integer exponents and 
low-lying excitations of different types, i.e. spin and charge with different 
velocities constituting so-called spin and charge separation, see [8-10]. 

A very curious situation arises in the context of the calculation of the par- 
tition function from the spectrum of the Hamiltonian. Despite the validity of 
the Bethe ansatz equations for all energy eigenvalues of the above mentioned 
models the direct evaluation of the partition function is rather difficult. In 
contrast to ideal quantum gases the eigenstates are not explicitly known, the 
Bethe ansatz equations provide just implicit descriptions that pose problems 
of their own kind. Yet, knowing the behaviour of quantum chains at finite 
temperature is important for many reasons. As a matter of fact, the strict 
ground-state is inaccessible due to the very fundamentals of thcrmodynam- 
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ics. Therefore the study of finite temperatures is relevant for theoretical as 
well as experimental reasons. 

The purpose of this review is to introduce to a unified treatment of the 
ground-state and the finite temperature properties of integrable quantum 
chains. First, in Sect. 2 we introduce the essential concepts of quantum in- 
tegrability. We show how to construct an infinite set of conserved currents, 
i.e. local operators commuting with the Hamiltonian. At the heart of these 
constructions is the embedding of the Hamiltonian into a family of commut- 
ing row-to-row transfer matrices of a certain classical model. In Sect. 3 we 
derive a lattice path integral representation for the partition function of a 
rather large class of integrable Hamiltonians at finite temperature. Here we 
also introduce a very efficient transfer matrix method based on the quantum 
transfer matrix. In Sect. 4 we review the algebraic Bethe ansatz for the semi- 
nal model of the partially anisotropic spin- 1/2 Hciscnbcrg chain (XXZ chain) 
related to the classical six- vertex model on a square lattice. The results of this 
section are the construction of eigenstates and eigenvalues of the Hamiltonian 
as well as the quantum transfer matrix of the system. In Sect. 5 we demon- 
strate how to transform the large number of coupled Bethe ansatz equations 
into a simple finite set of non-linear integral equations. These equations are 
studied numerically in Sect. 6 where explicit results for the temperature de- 
pendence of specific heat and susceptibility of the spin-1/2 XXZ chain are 
given. Finally, in Sect. 7 we resume the line of reasoning developed for the 
proof of intcgrability and study the thermal conductivity of the XXZ chain 
a topic that is of considerable current interest. 

2 Integrable exchange Hamiltonians 

We begin with the definition and discussion of a general class of quantum 
chains with nearest-neighbour interactions based on (graded) permutations. 
Consider a one-dimensional lattice with L sites and periodic boundary condi- 
tions imposed. A q-state spin variable on is assigned to each site i. We gener- 
ally consider the situation where each spin a has its own grading, i.e. statistics 
number e a = (— = ±1. A spin a with e a = +1 (e a = —1) is called bosonic 
(fermionic) . The Hamiltonian of the "permutation model" can be introduced 
as 

L 

ff = EVi (!) 

i=l 

with the (graded) permutation operator P»,j+i 

Pi.i+i l«i • ' ' cnai+i ■ ■ ■ a L ) = (-l)^^+i |ai • • • a i+1 ai ■ ■ ■ a L ) , (2) 

where ^ ai a i+1 is 1 if both a>i and a^+i are fermionic, and otherwise. 

Model (1) is shown to be exactly solvable on the basis of the Yang-Baxter 
equation. Many well-known exactly solvable models are of type (1), e.g. the 
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spin- 1/2 Heisenberg chain with q = 2 and e x = e 2 = +1, the free fermion 
model with q = 2 and t\ = —e 2 = +1, and the supersymmetric i-Jmodel 
with q = 3 and ei = — e 2 = e 3 = +1. If to of the e's are +1 and n(= q — m) 
arc — 1, for example ei = • • • = e m = +1, e m+i = • • • e q = — 1, we call the 
model (to, n)-permutation model or just (to, n)-model. 

Before sketching the proof of integrability of the general permutation 
model we look closer at two important special cases. The (2,0)-model (q = 
2, ei = 62 = +1) is the spin-1/2 Heisenberg chain with Hamiltonian 

L 

H = 2j2Si-S i+1 +L/2, (3) 

i=l 

in terms of SU(2) spin-1/2 operators S. 

The (2, l)-modcl (q = 3, ei = £2 = +1 and e 3 = —1) is the supersymmet- 
ric t — J model with Hamiltonian (ignoring a trivial shift) 

H = -tJ2n4,* c i+h* + c] +lttr c jia )V + J^2(S j S j+1 - n jnj+1 /4), (4) 

with standard fermionic creation and annihilation operators and c, pro- 
jector V = Y\j(l — nj^riji) ensuring that double occupancies of sites are for- 
bidden, and 2t — J (with normalization t = 1). The supersymmetric t — J 
model was shown to be integrable [11, 12] by the well-known Bethe ansatz 
[13, 14]. The ground-state and excitation spectrum were investigated [15] and 
critical exponents calculated by finite-size scaling and conformal field theory 
studies [16,17]. The thermodynamical properties were studied in [5] by use 
of the thermodynamical Bcthe ansatz (TBA) and in [18, 19] by use of the 
quantum transfer matrix (QTM). 

For proving the integrability of the quantum system, a classical counter- 
part is defined on a two-dimensional square lattice of L x N sites, where we 
impose periodic boundary conditions throughout this paper. We assume that 
variables taking values 1,2, •• • , q are assigned to the bonds of the lattice. 
Boltzmann weights are associated with local vertex configurations a, /?, /x 
and v and are denoted by R^, see Fig. 1. The classical counterpart to (1) is 
the Perk-Schultz model [20] with the following Boltzmann weights 

R p„( u , v ) = SavSfo + (u - v) ■ ■ S a0 S„ v , (5) 

where u and v are freely adjustable "interaction" parameters assigned to 
the entire horizontal and vertical lines intersecting in the particular vertex 
under consideration. These weights satisfy the Yang-Baxter equation which 
we depict graphically in Fig. 2. 

In the remainder of this section we want to indicate: (i) how the classical 
model is related to the quantum chain, and (ii) why the models are inte- 
grable. These issues are best discussed in terms of the transfer matrices of 
the classical model 
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R p v ( u > v ) = a 



Fig. 1. Graphical depiction of the fundamental i?-matrix (5). The Boltzmann 
weight assigned to each vertex configuration (configuration of spin variables around 
a vertex) corresponds to a matrix element of the matrix 7?. 




w 



Fig. 2. Depiction of the fundamental Yang-Baxter equation (YBE) for the ii-matrix 
with the following graphical rules: Each bond carries a spin variable (ranging from f 
to q), each vertex (via Fig. 1) corresponds to a local Boltzmann weight R depending 
on the local spin configuration. The algebraic term corresponding to each graph of 
the equation is obtained by multiplying the Boltzmann weights corresponding to the 
(three) vertices and summing over spin variables on closed (inner) bonds. The values 
of the spin variables on open (outer) bonds are fixed and the assignment of these 
values is identical for both sides. The statement of the graphical equation is that 
both sides of the equation evaluate to the same expression for arbitrary, however 
identical assignment of spins to the open bonds on either side. A consequence of 
the YBE is the commutativity of the transfer matrix with respect to the spectral 
parameter, see Sect. 3 and Fig. 6. 



L 

{a} i=l 

where we consider the spectral parameters Vi on the vertical bonds as fixed, 
i.e. independent of u. In dependence on the spectral parameter u the object 
T(u) represents a family of commuting matrices (the proof will be given in a 
slightly more general setting in the subsequent section) 

T(v)T(w) = T(w)T(v), for arbitrary v,w. (7) 

Commutativity holds especially in the case of vanishing spectral parameters 
Vi = on the vertical bonds, which case we refer to as the row-to-row transfer 
matrix. For this we have the additional simple limiting behaviour 
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T(0) = translation (shift) operator = c lP , (8) 
= Hamiltonian = H, (9) 



-f-lnT(u) 
du 



u=0 



also known as the Hamiltonian limit (P: momentum operator) [21,22]. 

Apparently, Hamiltonians obtained as members of commuting families of 
operators possess (infinitely) many conserved quantities. Any element of the 
family (or higher order derivative) commutes with H. For the case of isotropic 
SU (m, n)-symmetric systems we have presented the typical proof of integra- 
bility based on classical models satisfying the YBE. There are many more 
models satisfying the YBE with different or reduced symmetries. A famous 
example of a system with reduced symmetry is the partially anisotropic spin- 
1/2 Heisenberg chain (also known as XX Z chain). The Hamiltonian and 
i?-matrix corresponding to this will be given in the next section. 

Above, we have only shown how to incorporate the momentum opera- 
tor and the Hamiltonian into the family of commuting operators T(u). As 
indicated, lnT(u) is a generating function for conserved quantities 

^'-{ly^L- < io > 

that appear to be sums of local operators (with or without physical relevance). 
Quite generally, the current is related to the thermal current, whose 

conservation implies non-ballistic thermal transport, see Sect. 7. 

Logically, the presentation of diagonalization of the Hamiltonian or, more 
generally, of the row-to-row transfer matrix should follow directly after the 
above discussion of integrability. We will see, however, that another class 
of transfer matrices is worthwhile to be studied. This is a class of commut- 
ing staggered transfer matrices (quantum transfer matrices) occurring in the 
study of the thermodynamics of the quantum system at finite temperature 
(next section). We therefore postpone the discussion of the algebraic Bethe 
ansatz for both classes of transfer matrices to Sect. 4. 



3 Lattice path integral and quantum transfer matrix 

At the beginning of this section we want to comment on the various tech- 
niques developed for the study of thermodynamics of integrable systems at 
the example of the simplest one, namely the spin-1/2 Heisenberg chain. 

The thermodynamics of the Heisenberg chain was studied in [1, 2, 23] by 
an elaborate version of the method used in [24] . The macro-state for a given 
temperature T is described by a set of density functions formulated for the 
Bethe ansatz roots satisfying integral equations obtained from the Bethe 
ansatz equations (18). In terms of the density functions expressions for the 
energy and the entropy are derived. The minimization of the free energy 
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functional yields what is nowadays known as the Thermodynamical Bethe 
Ansatz (TBA). 

There are two "loose ends" in the sketched procedure. First, the descrip- 
tion of the spectrum of the Heisenberg model is built on the so-called "string 
hypothesis" according to which admissible Bethe ansatz patterns of roots are 
built from regular building blocks. This hypothesis was criticized a number 
of times and led to activities providing alternative approaches to the finite 
temperature properties. 

The second "loose end" within TBA concerns the definition of the en- 
tropy functional. In [1, 2, 23, 24] the entropy is obtained from a combinatorial 
evaluation of the number of micro-states compatible with a given set of den- 
sity functions of roots. As such it is a lower bound to the total number of 
micro-states falling into a certain energy interval. However, this procedure 
may be viewed as a kind of saddle point evaluation in the highly dimensional 
subspace of all configurations falling into the given energy interval. Hence, 
the result is correct in the thermodynamic limit and the "second loose end" 
can actually be tied up. Interestingly, the "second loose end" of the TBA ap- 
proach was motivation for a "direct" evaluation [25] of the partition function 
of integrable quantum chains. A straightforward (though involved) calcula- 
tion leads to the single non- linear integral equation of [26]. 

In this section we want to introduce the approach to thermodynamics of 
integrable quantum chains that we believe is the most efficient one, namely 
the "quantum transfer matrix" (QTM) approach. The central idea of this 
technique is a lattice path-integral formulation of the partition function of 
the Hamiltonian and the definition of a suitable transfer matrix [2 7-3 5]. 

In order to deal with the thermodynamics in the canonical ensemble we 
have to deal with exponentials of the Hamiltonian H . These operators are 
obtained from the row-to-row transfer matrix T{u) of the classical model in 
the Hamiltonian limit (small spectral parameter u) (9) 

T(u) = c ip + uH + ^ 2 \ (1) 

with P denoting the momentum operator. 

The main idea of the quantum transfer matrix (QTM) method at finite 
temperature is as simple as follows (for details the reader is referred to the 
papers [34, 35]). First, let us define a new set of vertex weights R by rotating 
R by 90 degrees as 

R;:(u,v) = Rf a (v,u), (2) 

see Fig. 3. We further introduce an "adjoint" transfer matrix T(u) as a prod- 
uct of R(— u, 0) [7, 34] with Hamiltonian limit 

T(u)^ C - ip+uH +°^ 2 \ (3) 

With these settings the partition function Zl of the quantum chain of length 
L at finite temperature T reads 
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— 0C|X 

R p v ( u ' v ) = a 



R p v (u,v) = a 



Fig. 3. Graphical depiction of the R and R matrices ((2) and (9)) associated with R. 
The Boltzmann weights assigned to the vertex configurations R and R correspond 
to the Boltzmann weights of the rotated fundamental vertices R. 



Z L = Tre- 0H = lim Z L , N (4) 

JV->oo 

where (3 = 1/T and Z^^ is defined by 

Z L , N := Tr [T(-r)T(-r)] N/ \ r := |. (5) 

The r.h.s.of this equation can be interpreted as the partition function of a 
staggered vertex model with alternating rows corresponding to the transfer 
matrices T(— r) and T(— r), see Fig. 4. We are free to evaluate the partition 
function of this classical model by adopting a different choice of transfer 
direction. A particularly useful choice is based on the transfer direction along 
the chain and on the corresponding transfer matrix T^™ which is defined 
for the columns of the lattice. The partition function of the quantum chain 
at temperature 1//3 is given by 

Z LtN = Tr (TQ™) L . (6) 

In the remainder of this paper we will refer to T QTM as the "quantum transfer 
matrix" of the quantum spin chain, because T Q ™ is the closest analogue to 
the transfer matrix of a classical spin chain. Due to this analogy the free 
energy / per lattice site is given just by the largest eigenvalue ^lmax of the 
QTM 

/ = -k B T lim log/l max . (7) 

N^oo 

Note that the eigenvalue depends on the argument r = (3/N which vanishes 
in the limit N — > oo requiring a sophisticated treatment. 

The QTM as defined above is actually a member of a commuting family 
of matrices T Q ™ defined by 

N/2 

[T Q ™] » = £ ft Ktr^T 1 («- r)Kl^l 2j (v, -r), (8) 

a j=l 

where we introduced yet another vertex weight J? as a rotation of R by -90 
degrees 
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N 



Fig. 4. Illustration of the two-dimensional classical model onto which the quantum 
chain at finite temperature is mapped. The square lattice has width L identical to 
the chain length, and height identical to the Trotter number N. The alternating rows 
of the lattice correspond to the transfer matrices T(—t) and T(— r), r = /3/N. The 
column-to-column transfer matrix T®™ (quantum transfer matrix) is of particular 
importance to the treatment of the thermodynamic limit. The arrows placed on 
the bonds indicate the type of local Boltzmann weights, i.e. R and /^-matrices 
alternating from row to row. (The arrows indicate the type of Boltzmann weight, 
they do not denote local dynamical degrees of freedom.) 



RZ(u,v) = R^(v,u). (9) 

Here we have introduced a spectral parameter v such that 7qtm(v) is a 
commuting family of matrices generated by v. A proof of this consists of 
two steps. First, we observe that R matrices and R matrices share the same 
intertwiner, i.e. the order of multiplication of two R matrices is interchanged 
by the same R matrix as in the case of the fundamental YBE, see Fig. 5. As R 
and R matrices share the same intertwiner, the transfer matrix obtained from 
products of arbitrary sequences of R and R with same spectral parameter (say 
v) on the continuous line constitutes a family of commuting matrices (spanned 
by v). The proof of this statement is graphically depicted in Fig. 6 for the 
case of the matrix T QTM (i>) (8). In other words, T QTM (w) is integrable which 
allows us to diagonalize Tqtm(w)- The final results, of course, are physically 
interesting just for v — as the physically meaningful QTM is identical to 
TQ™(0). 

The main difference to the transfer matrix treatment of classical spin 
chains is the infinite dimensionality of the space in which T Q ™ is living (for 
N — ► oo). In formulating (7) we have implicitly employed the interchange- 
ability of the two limits (L, N — ► oo) and the existence of a gap between 
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R(u,w)R(v,w) R(v,u)=R(v,u) R(v,w)R(u,w) 




R(u,w) R(u,v)R(w,v)=R(w,v)R(u,v) R(u,w) 

Fig. 5. Graphical derivation of the Yang-Baxter equation for R matrices. In the 
upper row, the fundamental YBE for R matrices is shown. The YBE for R matrices 
is obtained through rotation. Note that the intertwiner for R vertices is identical 
to the intertwiner for R vertices. 





Fig. 6. Railroad proof for the commutativity of transfer matrices: The intertwiner 
is pulled through from right to left by a successive application of the YBE. 



the largest and the next-largest eigenvalues of T Q ™ for finite temperature 
[36,37]. 
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The next-leading eigenvalues give the exponential correlation lengths £ of 
the equal time correlators at finite temperature 



- = lim In 



(10) 



Lastly we want to comment on the study of thermodynamics of the quan- 
tum chain in the presence of an external magnetic field h coupling to the 
spin S = Y^j=i Sji where Sj denotes a certain component of the jth spin, for 
instance Sj. This, of course, changes (5) only trivially 

Z L , N := Tr { [T{-r)T(-r)] m ■ e? hs } . (11) 

On the lattice, the equivalent two-dimensional model is modified in a simple 
way by a horizontal seam. Each vertical bond of this seam carries an individ- 
ual Boltzmann weight e ±/3/l / 2 if Sj = ±1/2 which indeed describes the action 
of the operator 

L 

( J >hS = Y[<J >hS i. (12) 
j'=i 

Consequently, the QTM is modified by an h dependent boundary condition. 
It is essential that these modifications can still be treated exactly as the 
additional operators acting on the bonds belong to the group symmetries of 
the model. 

We like to close this section with some notes on the relation of the two 
apparently different approaches, the combinatorial TBA and the operator- 
based QTM. In fact, these methods are not at all independent! In the latter 
approach there are several quite different ways of analysis of the eigenvalues 
of the QTM. In the standard (and most economical) way, see below, a set 
of just two coupled non-linear integral equations (NLIE) is derived [34,35]. 
Alternatively, an approach based on the "fusion hierarchy" leads to a set of 
(generically) infinitely many NLIEs [34, 38] that are identical to the TBA 
equations though completely different reasoning has been applied! 

Very recently [26] , yet another formulation of the thermodynamics of the 
Heisenberg chain has been developed. At the heart of this formulation is just 
one NLIE with a structure very different from that of the two sets of NLIEs 
discussed so far. Nevertheless, this new equation has been derived from the 
"old" NLIEs [26,39] and is certainly an equivalent formulation. In the first 
applications of the new NLIE, numerical calculations of the free energy have 
been performed with excellent agreement with the older TBA and QTM 
results. Also, analytical high temperature expansions up to order 100 have 
been carried out on this basis. 
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4 Bethe ansatz equations for the spin-1/2 XX Z chain 

In the following we consider the anisotropic Heiscnberg chain (slightly gen- 
eralizing (3)) with Hamiltonian Hl 

H = 2"£[SfSf +1 + Sp y 1+1 + AS*S* +i } (1) 

with periodic boundary conditions on a chain of length L. Apparently, for 
A = +1 the system specializes to the isotropic antifcrromagnctic Heiscnberg 
chain, for A = — 1 (and applying a simple unitary transformation) the system 
reduces to the isotropic ferromagnetic case. 

The classical counterpart of the XX Z chain is the six- vertex model. For 
our purposes the following parameterization of the Boltzmann weights is 
useful 

, x , j . , sin(7w/2) . , sin7 , . 

a{w) = 1, b(w) = , - K '' ' - , c(w) = — jl -. (2) 

K ' K ' sin(7w/2 + 7 )' v ' sin( 7 w/2 + 7) v ; 

All relations between Hamiltonian of the quantum system, row-to-row trans- 



a(u-v) 



b(u-v) 



c(u-v) 



Fig. 7. Allowed vertices of the six-vertex models with corresponding Boltzmann 
weights a, b, c. The remaining vertices are obtained by an exchange of bond-states 
1«2. 



for matrix of the classical model and the QTM as introduced in the previous 
section are also valid in the present case if A and 7 are related by A — cos 7. 
The only modification concerns (9) as here the relation of Hamiltonian (1) and 
the logarithmic derivative of the row-to-row transfer matrix for the weights 
(2) aquires a normalization factor 

F = 2^A lnT ( u )| (3 ) 
7 du x ;]u=0 K 1 

Monodromy matrix 

Our aim is to diagonalize the row-to-row transfer matrix and the QTM by 
means of the algebraic Bethe ansatz. We first review some notation and basic 
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properties of i?-matrices (as collections of local Boltzmann weights) and the 
so-called L-matrix. The elements of the L-matrix at site j are operators acting 
in the local Hilbert space hj (for the six- vertex model isomorphic to ~ C 2 ). 
The L-matrix' element in row a and column j3 is given in terms of the R 
matrix 



where is a matrix with only non-vanishing entry 1 in row \i and column 
v. This reads explicitly for the six-vertex model 



with a, b, c given in (2). 1 The so-called monodromy-matrix is defined as a 
product of all L-matrices on consecutive sites 



In essence, the (a, /3)-element of the monodromy matrix is the transfer ma- 
trix of a system with fixed boundary spins a on the left end and j3 on the 
right end of a row. 

Algebraic Bethe ansatz 

The procedure of diagonalization can be decribed as follows: 

First we search for a pseudo- vacuum state (reference state) \Q) that is 
a simple eigenstate of the operator-valued diagonal entries A and D of 
the monodromy matrix T (and hence an eigenstate of the transfer matrix 
T = A + D). The lower-left entry C of the monodromy matrix applied to 
\fi) yields zero, the upper-right entry B yields new non- vanishing states. 
Hence C and B play the role of annihilation and creation operators. 

From the Yang-Baxter equations a quadratic algebra of commutation re- 
lations for the entries of the monodromy matrix (notably for A, D, and 
B) is derived. By use of these relations algebraic expressions of the eigen- 
states and rather explicit expressions for the eigenvalues are derived. 

It turns out that the row-to-row transfer matrix and the QTM can be 
treated in parallel as the quadratic algebra of the entries of the monodromy 
matrix is identical for both cases. The only difference lies in the different 
reference states and the different "vacuum expectation values" . We therefore 
focus our presentation on the slightly simpler case of the row-to-row transfer 

1 As a reminder: a + — a x + ia v and er~ = a x — ia v , and states |+), | — ) correspond 



L j «{w)=Rj v {w)e j », 



(4) 







to |1>, |2>. 
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matrix and describe the necessary modifications for the QTM in (20,21). 
Reference state 

Provided we find local states \uii) that are eigenstates of the diagonal entries 
of Li and are annihilated by the left-lower entry, then | Q) may be taken as 
the product of these local states. A glance to (5) shows that this is possible 
just with \u)i) = \2) i 

N 

l«)=®|2) i . (7) 

i 

The monodromy matrix T applied to | Q) yields an upper triangular 2x2 
matrix of states 

(a N \Q) B\Q) 
J W-{ b N \Q) 

or explicitly 

A\Q) =a N \f2), D\f2)=b N \[2), T \f2) = (a N + b N ) \f2) . (9) 
Therefore, J?) is an eigenstate of T. 

Quadratic algebra of operators A, B, and D 

We intend to use the operator B as creation operator for excitations, i.e. we 
demand that the new state \Q\(v)) := B(v) \Q) ("one-particle state") be an 
eigenstate of T(u) = A(u) + D{u). What we need to know is the operator 
algebra for interchanging B(v) with A(u) and D(u). This algebra can be 
obtained from the YBE, for a graphical representation see Fig. 8. By fixing 




u 




v — u = v — u 




T(u)®T (v)R(v -u) = R(v - u)T (v) <g> T (u) 



Fig. 8. Graphical depiction of the YBE. 



the exterior spins on the horizontal bonds we obtain all the commutators we 
need for interchanging the operators. We begin with the relation of any two 
£>-operators illustrated in Fig. 9 actually implying commutation 



[B(u),B(v)}=0. 



(10) 
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- ^y < c 



B(u)B(v)a(v -u) = a(v - u)B(v)B{u) 



Fig. 9. Quadratic relation of B-operators obtained from fixing the spins on the left 
and right open bonds in Fig. 8. Note that the summation over the spin variables 
on the closed bonds between R and T matrices reduces to just one non-vanishing 
term on each side of the equation. 



Now we look at products of A and B and again use a graphical represen- 
tation shown in Fig. 10. We algebraically find 



A(u)B(v) = ^ U } b(v)A(u) - fj" U \ b(u)A(v). 



b(v — u) 



b(v — u) 



(11) 



2/1 





1/2 " 1 

B(u)A(v)c(v -u) + A{u)B(v)b(v - u) = a(v - u)B(v)A(u) 



Fig. 10. Bilinear relation of B and A-operators obtained from fixing the spins 
on the left and right open bonds in Fig. 8. Note that the summation over the 
spin variables on the closed bonds between R and T matrices reduces to two non- 
vanishing terms on the left side of the equation and to one non- vanishing term on 
the right side of the equation. 



In a similar same way, illustrated in Fig. 11, wc obtain (after interchanging 

u «-> v ) 

D(u)B{v) = a J; U ~ V l B{v)D(u) - C ^~ V \ b{u)D{v). (12) 



b(u — v) 



b(u — v) 



Now we have available all necessary relations for evaluating the applica- 
tion of T(u) to By use of the commutation relations we get 
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V 











B(u)D(v)a(v -u)= c(v - u)B(v)D{u) + b(v - u)D(v)B(u) 

Fig. 11. Bilinear relation of B and D-operators obtained in analogy to the above 
cases. 



T{u)\Q,{v)) 



, , a(v — u) a(u - v) 



a(v) 



b(v — u) 
c(v — u) 
b(v — u) 



b(u — v) 
c(u — v) 
b(u — v) 



B(u) \Q) 



(13) 



where we have used the abbreviations a = a N , f3 = b N . The terms in the 
first (second) line of the r.h.s.of (13) are called "wanted terms" ("unwanted 
terms" ) and arise due to the first (second) terms on the r.h.s.of (11,12). In 
order that \£2i(v)} be an eigenstate of the transfer matrix the second term in 
(13) has to vanish and the first one to give the eigenvalue 2 



a(v) c{u — v)b(v — u) 



(3(v) " 

A{u) = a(u) 



c(v — u)b(u — v 
a(v — u) 



b(v — u) 



= 1, 
a(u — v) 



b(u — v) 



(14) 
(15) 



We can generalize this to any n-particle state. The argument is quite the 
same as for just one excitation. We look at the following state 



|fl(«i)> = tlB(vi)\n) 



(16) 



where the numbers Vi will be referred to as Bethe ansatz roots. Demanding 
that this state be an eigenstate of T(u) we get after successive application of 
the commutation rules the following set of equations. From the two "wanted 
terms" the eigenvalue is read off 



A{u) = a{u) J] 



(17) 



and vanishing of the "unwanted terms" yields 

2 For the sake of transparency of the algebraic structure we do not use the fact 
that by our convention the Boltzmann weight a is identical to 1. 
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a{vi) 



n 



K v j - v i) 
b{vi-VjY 



for i = 1, 



(18) 



The last constraints are nothing but the famous Bethe ansatz equations. 
We like to note that we would have obtained the same set of equations by 
demanding that the function on the r.h.s.of (17) be analytic in the whole 
complex plane 

ResA(u = Vi)=0 Mi. (19) 

In the case of a general transfer matrix T (v; {u^}; {u^}j defined as a 

product of Ni many vertices of type R and N 2 many vertices of type R with 
spectral parameter v in the auxiliary space (corresponding to the continuous 
line) and in the quantum spaces the eigenvalue expression is very 

similar to (17) and (18). The only change is the replacement 



N 2 



a(u) 



(u — u 



(2) 



; +^/ 2 n o(t 

fe=l 

iVi 



u), 



k=l 
N 2 



(20) 



(21) 



k=l 



k=l 



where we have also introduced exp(±/3/i/2) factors arising from twisted 
boundary conditions for the transfer matrix as realized in the case of the 
quantum transfer matrix for a system in a magnetic field h. 

We may simplify the following presentation if we perform a rotation by 
7r/2 in the complex plane, i.e. we introduce a function A by X(v) = A(iv) 



and for convenience we replace Vj 



The eigenvalue expression for any 



eigenvalue X(v) of the transfer matrix T reads (see also [22] 

a(ivj — if) ( ( . ^ yj a(iv 



X(v) = a(iv) [] 



\ Kivj - iv) 



(22) 



Of particular importance is the case of the quantum transfer matrix with 



,(2) 



Ni = N 2 = N/2 and all 
normalization factor in (3). 
Introducing the definition 

s(v) :— sinh(7v/2), 

we obtain for the functions b, a, [3 

s(v) _ 
i)' 

x -| N/2 



with 



2^4, note the 

7 N 1 



(23) 



b(iv) = 



s(v - 2i)' 



a(iv) = e^/ 2 



s(v 



s(v 



2i) 



0(iv) 



-0h/2 



s(v + it) 
s(v + it — 2i) 



N/2 



(2.4) 
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From (22) we obtain for the function X(v) 

x(v ) = AiM + W (25) 

{> [s(v-i(2~T))s{v + i(2-T))] N /2' V > 
where the terms Xi^(v) are 



q(v) 

A 2 («):=e-^/ 2 ^ + i)fc^), (26) 

and <j>{v) is simply 



<j>{v) := [s(v - i(l - t))s(v + i(l - r))]^ 2 . (27) 

The function q(v) is defined in terms of the yet to be determined Bethe ansatz 
roots vj 

g(«):=n*(«-"i)- (28) 
i 

Note that we are mostly interested in A which is obtained from X(v) simply by 
setting v = 0. Nevertheless, we are led to the study of the full ^-dependence 
since the condition fixing the values of vj is the analyticity of Xi(v) + \2(v) 
in the complex plane. This yields 

a( Vj ) = -1, (29) 

where the function a(v) (not to be confused with the Boltzmann weight a(w) 
above) is defined by 

a{v) := ^4 = e ^ y-j)g^+2i) (30) 

v ; A 2 (u) ^(w + i)q(v - 2i) v 7 

Algebraically, we are dealing with a set of coupled non-linear equations sim- 
ilar to those occurring in the study of the eigenvalues of the Hamiltonian 
[22]. Analytically, there is a profound difference as here in (30) the ratio of 
(^-functions possesses zeros and poles converging to the real axis in the limit 
N — > oo. As a consequence, the distribution of Bethe ansatz roots is discrete 
and shows an accumulation point at the origin, cf. Fig. 12. Hence the treat- 
ment of the problem by means of linear integral equations for continuous 
density functions [40] is not possible in contrast to the Hamiltonian case. 



5 Manipulation of the Bethe ansatz equations 

The eigenvalue expression (25) under the subsidiary condition (29) has to 
be evaluated in the limit N — > oo. This limit is difficult to take as an in- 
creasing number A^/2 of Bethe ansatz roots Vj has to be determined. In the 
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Fig. 12. Sketch of the distribution of Bethe ansatz roots Vj for finite N. Note that 
the distribution remains discrete in the limit of N — > oo for which the origin turns 
into an accumulation point. 



Hamiltonian case, i.e. for (20) with N\ = N and N 2 — (or vice versa), 
the distribution of roots is continuous and the Bethe ansatz equations (18) 
can be reduced to linear integral equations [22, 40]. For the general case, this 
is no longer possible. In this case, notably for the QTM, the distribution of 
the roots is discrete and the standard approach based on continuous density 
functions is not possible. From now on we explicitly discuss the QTM case, 
i.e. (20) with Ni = N 2 = N/2 yielding all information on the free energy 
at arbitrary temperature T including the limit T = 0. (With some modifi- 
cations, the final non-linear integral equations also apply in the case of the 
row-to-row transfer matrix and the Hamiltonian. These results, however, are 
not of prime interest to this review.) 



5.1 Derivation of non-linear integral equations 

The main idea of our treatment is the derivation of a set of integral equations 
for the function a(v). This function possesses zeros and poles related to the 
Bethe ansatz roots Vj, see Fig. 13. Next we define the associated auxiliary 
function A(v) by 

A(v) = l + a(v). (1) 

The poles of A(v) are identical to those of a(v). However, the set of zeros is 
different. From (29) we find that the Bethe ansatz roots are zeros of A(v) 
(depicted by open circles in Fig. 14). There are additional zeros farther away 
from the real axis with imaginary parts close to ±2. For the sake of complete- 
ness, these zeros are depicted in Fig. 14 (open squares), but for a while they 
are not essential to our reasoning. Next we are going to formulate a linear 
integral expression for the function loga(w) in terms of log A(v). To this end 
we consider the function 

1 f d 

f( v ) : =7r^ / -r log s(v - w) log A(w)dw (2) 
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Fig. 13. Distribution of zeros (o) and poles (x) of the auxiliary function a(v). All 
zeros and poles Vj =p 2i are of first order, the zeros and poles at ±(2i — ir), ±ir arc 
of order N/2. 
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Fig. 14. Distribution of zeros and poles of the auxiliary function A(v) = 1 + a(v). 
Note that the positions of zeros (o) and poles (x) are directly related to those 
occurring in the function a(v). There are additional zeros (□) above and below the 
real axis. The closed contour C by definition surrounds the real axis and the zeros 
(o) as well as the pole at — ir. 



defined by an integral with closed contour C surrounding the real axis, the 
parameters vj and the point — ir in anticlockwise manner, see Fig. 14. Note 
that the number of zeros of A(v) surrounded by this contour is N/2 and 
hence is identical to the order of the pole at — ir. Therefore the integrand 
log A(w) does not show any non-zero winding number on the contour and 
consequentially the integral is well-defined. By use of standard theorems we 
see that the function f(v) is analytic in the complex plane away from the 
real axis (and axes with imaginary parts tiples of 2-77/7). Next we perform an 
integration by parts and apply Cauchy's theorem 
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f(v) = J2 l0 Ss(v - vj) - N/2 \ogs(v + ir) = log g( "L /2 , (3) 

j S(V + 1T) ' 

Thanks to (2) and (3) we have a linear integral representation of log q(v) in 
terms of log A(v) . Because of (30) the function log a(v) is a linear combination 
of log q and explicitly known functions leading to 

log a(v) =Ph + log ( f ""if + ^ i + iT n ^ + f(v + 2i) -f(v- 2i). (4) 



v s(w + it)s(v + 2i — ir) / 

From now on we use a shorthand notation for the logarithmic derivative of 

s(v) 

d 7 7 

d(v) := — log s(v) = - coth -v. (5) 

By use of the definition of the integration kernel k 

1 d s(u — 2i) 1 7 sin 27 



x( u ) : = — — log -) -( = — [d(u - 2i) -d(u + 2i)l = „ 

v ; 2m du s(w + 2i) 2m l y ' K n 2n cosh -yu - cos 2-y 

(6) 

we may write (4) as 

N i ( s(v-iT)s(v + 2i + iT)\ f . , n A , 

log a(v) = $h + — log ( ) - / log i w dw. 

2 V s ( w + it)s(u + 2i - it)/ J £ 

(7) 

This expression for a(v) is remarkable as it is a non-linear integral equation 
(NLIE) of convolution type. It is valid for any value of the Trotter number 
N which only enters in the driving (first) term on the r.h.s.of (7). This term 
shows a well-defined limiting behaviour for TV — > 00 



N f s(v -it)s(v + 2i + it) , 
— log — — — — > iNt 

2 6 \s(v + ir)s(v + 2i - it) ' 



d , , d , , . 

— log s(u + 2i) — log s(w) 

aw dv 



= i(32^1[d(v + 2i)-d(v)}, (8) 

7 

leading to a well-defined NLIE for a(v) even in the limit N — > 00 

\oga(v) — ph + (3e (v + i) — J >c(v — w) log A(w)dw, (9) 

where eo is defined by 

e («) = + i) - d(« - i)] = 2 Sm 7 . (10) 

cosh yv - cos 7 

This NLIE allows for a numerical (and in some limiting cases also analytical) 
calculation of the function a(v) on the axes = ±1. About the historical 
development we like to note that NLIEs very similar to (9) were derived 
for the row-to-row transfer matrix in [41,42]. These equations were then 
generalized to the related cases of staggered transfer matrices (QTMs) of the 
Heisenberg and RSOS chains [34,35] and the sine-Gordon model [43]. 
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5.2 Integral expressions for the eigenvalue 

In (7) and (9) we have found integral equations determining the function a 
for finite and infinite Trotter number N, respectively. The remaining problem 
is the derivation of an expression for the eigenvalue A in terms of a or A. 

From (25) we see that X(v) is an analytic function with periodicity in 
imaginary direction (period 27ri/7) and exponential asymptotics along the 
real axis. Hence, up to a multiplicative constant C, we can write 

Xi(v)+X 2 (v) = Cl[s(v-w l ), (11) 

i 

where w;,/ = 1, N, are the zeros of X(v) — Ai(w)+A 2 (w) which are solutions 
to a(v) = Xi(v)/X 2 (v) = — 1, i.e. zeros of A(v) = l + a(v) that do not coincide 
with Bethe ansatz (BA) roots! These zeros are so-called hole-type solutions 
to the BA equations. The holes are located in the complex plane close to the 
axes with imaginary parts ±2, see zeros in Fig. 14 depicted by □. Thanks to 
Cauchy's theorem we find for v in the neighbourhood of the real axis 



/ d(v-w- 2i)[log A(w)]'dw = ^d(v- vj - 2i) - ^-d{v + ir - 2i) 
71-1 c 3 

(12) 

as the only singularities of the integrand surrounded by the contour £ are 
the simple zeros vj and the pole — ir of order N/2 of the function A. Also, 
for v in the neighbourhood of the real axis we obtain 



[ d(v - w) [log A(w)]'dw = d(v - vj - 2i) - ^ d(v 
71-1 c j i 



wi) 



3 

N 

+ -d(v + 2i-iT), (13) 

where the evaluation of the integral has been done by use of the singularities 
outside of the contour £ and use of the period 2iri/j of the integrand. The 
old contour £ is replaced by a contour £ such that the upper (lower) part of 
£ is the lower part of £ (the upper part of £ shifted by — 27ri/7), and reversed 
orientation. The surrounded singularities are the simple poles Vj + 2i — 27ri/7, 
the zeros wi with or without shift —2tti/j, and the pole ir — 2i of order N/2 
of the function A. 

Next, we take the difference of (12) and (13), perform an integration by 
parts with respect to w, and finally integrate with respect to v 

— / [d(v — w) — d(v — w — 2i)] log A(w)dw 
2tt\ J c 

_ [.(.-.(» -r)M. + i(2-r))]^ + const 

Hi s ( v - w o 

Combining (11), (14) and (25) we find 
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logA(u) = ~(3h/2 - — [ [d(v -w)-d(v-w- 2i)] log A(w)dw, (15) 
2ni J c 

where the constant was determined from the asymptotic behaviour for v — > oo 
and use of A(oo) = exp(/3/i/2) + exp(— (3h/2) and A(oo) = 1 + exp(/3/i). 

These formulas, (15) and (9), are the basis of an efficient analytical and 
numerical treatment of the thermodynamics of the Heisenberg chain. There 
are, however, variants of these integral equations that are somewhat more 
convenient for the analysis, especially for magnetic fields close to 0. The 
alternative integral expression for A reads [34, 35] 

r°° — i 

In A = -0eo + / K(x) \n[Vl(x)%(x)}dx, K(x) = - — — - , (16) 

where 2l(x) and 2l(x) are complex-valued functions with integration paths 
along the real axis. These functions are determined from the following set of 
non-linear integral equations 



lna(x) = -/3 '- — — + — + K*]n%L(x) -k* In 21 (x + 2i), 

7 cosh I a; 2(n - 7) 

sin t 7T irSh 

lna(x) = -/3 '- — — + k * lnQllx) -K*ln%.(x- 2i) 

y ' ' 7 cosh fx 2(tt-7) K ' V ' 



(17) 



(18) 

2l(x) = 1 + a(x), 2l(x) = 1 + a(x). (19) 

The symbol * denotes the convolution / * g(x) = f(x — y)g(y)dy and 
the function k(x) is defined by 



1 roc sinh ( — — 2 ) k 

< x ) = TT" / V7 -, J ^e ikx dk. (20) 

2?r 7_oo 2 cosh k sinh ( Z - 1 J k 



Note that the integrals in (17) and (18) are well-defined with integration 
paths just below and above the real axis. 

The above equations are obtained from (9) by a partial "particle-hole" 
transformation of the function a(v) only on the axis 3(v) = — 1. Replacing 
log A = log A + log a (where a — 1/a, A = 1 + a) on the lower part of C in 
(9) leads to an equation involving convolution type integrals with log A, log A 
and log a. This equation can be resolved explicitly for log a by straightforward 
calculations in "momentum space". Finally, a(x) := a(x + i) and a(x) := 
a(x — i). 



6 Numerical results for thermodynamical quantities 



By numerical integration and iteration the integral equation (9) can be solved 
on the axes 3(w) = ±1 defining functions a ± (x) := a(x±i). Alternatively and 
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particularly convenient for the case of vanishing magnetic field h, equations 
(17) and (18) can be used for the functions a and a. Choosing appropriate ini- 
tial functions the series a k with k = 0, 1, 2, . . . converges rapidly. In practice 
only a few steps are necessary to reach high-precision results. Moreover, using 
the well-known Fast Fourier Transform algorithm we can compute the convo- 
lutions very efficiently. In fact, some of the convolutions in (9) or (17,18) are 
delicate to be evaluated in "real space" , because of the appearance of a pole of 
the kernel just at the integration contour. These integrals are automatically 
handled correctly in "momentum space" . 

In order to calculate derivatives of the thermodynamical potential with 
respect to temperature T and magnetic field h one can avoid numerical dif- 
ferentiations by utilizing similar integral equations guaranteeing the same 
numerical accuracy as for the free energy. The idea is as follows. Consider 
the function 

1 da ci 

lap, 



d d 
lap:= log a with — log(l + a) 



I + ad/3 l + i 



we directly obtain from (9) a linear integral equation for lap if we regard 




Fig. 15. Specific heat c(T) data versus temperature T for the spin-1/2 XX Z chain 
with repulsive interaction < A < 1. 



the function a as given. Once the integral equation (9) is solved for a, the 
integral equation for lap associated with (9) can be solved. In this manner, we 
may continue to any order of derivatives with respect to T (and h). However, 
in practice only the first and second orders matter. Here we restrict our 
treatment to the specific heat c(T) and the magnetic susceptibility x(T) 
(derivatives of second order with respect to T and h), see Figs. 15 and 16. 
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Fig. 16. Susceptibility x(^) data versus temperature T for the spin-1/2 XX Z 
chain with repulsive interaction < A < 1. 



Note the characteristic behaviour of c(T) and x(T) at low temperatures. 
The linear behaviour of c(T) and the finite ground-state limit of x(T) are 
manifestations of the linear energy-momentum dispersion of the low-lying ex- 
citations (spinons) of the isotropic antifcrromagnctic Hcisenbcrg chain. Also, 
with increase of the repulsion A the location of the finite temperature max- 
imum in c(T) shifts to higher temperatures and the values of x(T) drop. In 
the high temperature limit the asymptotics of c(T) and x(T) are 1/r 2 and 
1/T. This and the existence of the finite temperature maximum are a con- 
sequence of the finite dimensional local degree of freedom, i.e. the spin per 
lattice site. 

For A = 1 note that x(T) approaches the ground-state limit x(0) m a 
singular manner, see also Fig. 17. The numerical data at extremely low tem- 
peratures provide evidence of logarithmic correction terms, see also [44] and 
later lattice studies [45,46] confirming the field theoretical treatment by [47]. 
These logarithmic terms are responsible for the infinite slope of x(^) at T = 
despite the finite ground-state value x(0) = l/V 2 . Precursors of such strong 
slopes have been seen in experiments down to relatively low temperatures, 
see e.g. [48]. Unfortunately, most quasi one-dimensional quantum spin sys- 
tems undergo a phase transition at sufficiently low temperatures driven by 
residual higher dimensional interactions. Hence the onset of quantum critical 
phenomena of the Heisenberg chain at T — becomes visible, but cannot be 
identified beyond all doubts. 
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Fig. 17. Magnetic susceptibility x at low temperature T for the isotropic spin-1/2 
XXX chain. In the inset x{T) is shown on a larger temperature scale. 



7 Thermal transport 

The simplest approach to the investigation of transport properties is based on 
linear response theory leading to the Kubo formulas [49, 50] relating conduc- 
tivities to dynamical correlation functions of local current operators. Hence, 
the calculation of transport properties is more difficult than the computa- 
tion of the free energy. In fact, the calculation of correlation functions is a 
most difficult issue and exact results are rare, even for integrable systems. In 
some situations, however, the explicit computation of correlation functions 
can be avoided. For the spin transport of the Heisenberg chain (= electrical 
transport of the model in the particle representation) the Drude weight of 
the dynamical conductivity can be cast into a form involving only spectral 
properties without explicit recourse to matrix elements. At the time of writ- 
ing, these expressions are known how to be evaluated for zero temperature, 
the case of non-zero temperature is still controversial. 

A fortuitous case is the thermal transport as the thermal current Je is 
one of the conserved currents (10). The first three conserved currents 

(n = 0, 1, 2) are related to the momentum operator, the Hamiltonian and the 
thermal current via 
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P = -ij (0) (1) 
7 2 

J E = i(2^ 2 J^ + iL. (2) 

Note that the spin current (electrical current) is not conserved as it is not 
contained in the sequence of conserved currents ! 

Let us first motivate that is indeed the thermal current Je of the 
system. To this end we impose the continuity equation relating the time 
derivative of the local Hamiltonian (interaction) and the divergence of the 
current: h = — divj E . The time evolution of the local Hamiltonian hk : k+i is 
obtained from the commutator with the total Hamiltonian and the divergence 
of the local current on the lattice is given by a difference expression 

^giW = i[Ht hkk+l{t)] = _ {ifeE+l(t) _ fi m (3 ) 

Apparently the last equation is satisfied with a local thermal current operator 
jf defined by 

jk = i[hk-ik,hkk+i]- (4) 

In fact, up to a trivial scale and shift, the conserved current is identical 
to the r.h.s.of the upper equation, see also [51-56]. 

The Kubo formulas [49, 50] are obtained within linear response theory 
and yield the (thermal) conductivity k relating the (thermal) current Je to 
the (temperature) gradient VT 

Je - kVT, (5) 

where 

poo pP 

k(l>)=0 dte-' lu,t dT(M-t-iT)J E ). (6) 
Jo Jo 

As the total thermal current operator Je commutes with the Hamiltonian H 
of the XX Z chain we find 



1 

i(w — ie) ' 



«M = (e-0+), (7) 



with 5Rk(w) = k5(gj) where 

k = irp 2 (Ji). (8) 

As a consequence, the thermal conductivity at zero frequency is infinite! This 
is only natural, as the conserved current cannot decay in time. However, the 
weight of the zero-frequency peak is some finite and non-trivial temperature 
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dependent quantity to be calculated from the second moment of the thermal 
current. 

Quite generally, the expectation values of conserved quantities may be 
calculated by use of a suitable generating function 



Z = Tr exp(-/3H + XJe), 



(9) 



from which we find the expectation values by derivatives with respect to A 
at A = 



A=0 



■Je) =0. (A) lu2 



A=0 



(Ji)-(J E ) 2 = (Ji), (io) 



where we have used that the expectation value of the thermal current in 
thermodynamical equilibrium is zero. 

Instead of Z we will find it slightly more convenient to work with a par- 
tition function 

Z = Trexp(-A 1 j( 1 )-A„j(")), 



notably n = 2. With view to (2) we choose 

sin 7 



Ai = 2 



7 



A 2 = -iA 2 



sin 7 



and obtain the desired expectation values from Z 



A=0 



(11) 



(12) 



(13) 



We can deal with Z rather easily. Consider the trace of a product of TV 
row-to- row transfer matrices T(uj) with some spectral parameters Uj close 
to zero, but still to be specified, and the Nth power of the inverse of T(0) 



Z N = Tr [T(«i) T(u N ) ■ T(0)- N ] 

= Trexp[^[lnT( % )-lnT(0)] I . 



(14) 



Now it is a standard exercise in arithmetic to devise a sequence of N numbers 



u\,...,Un (actually Uj = such that 

Ihx l £[/(«i)-/(0)] = -A 1 |-/(«) 



— \ r . 



u=0 



0_ 

du 



/(«) 



u=0 



(15) 



We only need the existence of such a sequence of numbers, the precise values 
are actually not important. In the limit N — > oo we note 
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lim Z N = Z. (16) 



We can proceed along the established path of the quantum transfer matrix 
(QTM) formalism presented above and derive the partition function Zn in 
the thermodynamic limit L — > oo 

lim Z]I L = A, (17) 

L — ^oo 

where A is the largest eigenvalue of the QTM. The integral expression for A 
reads 



/oo 
K(x) ln[%(x)%(x)]dx, 
-oo 



(18) 



with K{x) already defined in (16) and some function e(x) given in [34,35]. 
In the limit N — ► oo the first term on the r.h.s.of the last equation turns into 

Ft / ft \ n 

Jim 5>( ttJ ) - e( 0)] = -Ai^eOO^ - A„ (-) e { u)\^ (19) 



a rather irrelevant term as it is linear in Ai and A„, and therefore the second 
derivatives with respect to Ai and A„ vanish. The functions 2l(x) and 2l(x) 
are determined from the following set of non-linear integral equations 

lna(x) = ^[e (a; - iuj) - £o(0)] + k * ln2l(x) — k * ln2l(x + 2i), 

3 

lna(x) = ^[e (x — iuj) — So(0)] + k * ln2l(x) — k * ln2l(x — 2i), 
j 

2l(x) = 1 + a(x), 2l(x) = 1 + a(x). (20) 

with a function So(x) given in terms of hyperbolic functions [34,35]. Again, 
the summations in (20) can be simplified in the limit N — > oo 

' » .. > » .. ' 



— :si(x) =:e n (x) 

(21) 

where the first function can be found in [34, 35] and is simply 



and hence the second function is 



d\ n— 1 



£n[x) = I ' (23) 
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We like to note that the structure of the driving term (21) appearing in 
the NLIE (20) reflects the structure of the generalized Hamiltonian in the 
exponent on the r.h.s.of (11). We could have given an alternative derivation 
of the NLIE along the lines of the thermodynamic Bethe ansatz (TBA). 
In such an approach the driving term is typically the one-particle energy 
corresponding to the generalized Hamiltonian. Hence it has contributions 
due to the first as well as the nth logarithmic derivative of the row-to-row 
transfer matrix, i.e. the terms e\ and s n . 

In Fig. 18 we show k(T) for various anisotropy parameters 7. Note that 
k(T) has linear T dependence at low temperatures. At high temperatures 
k(T) behaves like 1/T 2 . 
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Fig. 18. Illustration of numerical results for the thermal conductivity R as a 
function of temperature T for various anisotropy parameters A = cos 7 with 
7 = 0, tt/6, 7r/5, tt/4, tt/3, tt/2. 



8 Summary 

We have reviewed the treatment of intcgrablc quantum systems at zero and 
finite temperature based on a lattice path integral formulation and the def- 
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inition of the so-called quantum transfer matrix (QTM). In detail, the ther- 
modynamical properties of the partially anisotropic Heisenberg chain were 
discussed. As we hope, a transparent analysis of the eigenvalue problem of 
the row-to-row and quantum transfer matrices has been given, resulting in 
a set of non- linear integral equations (NLIE). From a numerical solution of 
these NLIEs at arbitrary temperature the specific heat and magnetic suscep- 
tibility data were obtained. Also the Drude weight of the thermal current 
was calculated. Although we have given the explicit results only for the case 
of vanishing magnetic field (h = 0), the derivation of the free energy is not 
restricted to this case. In fact, the NLIEs given above are valid even in the 
case of non- vanishing fields h. However, the above treatment of the thermal 
conductivity is limited to the case h = 0. The generalization to arbitrary 
fields is an open problem as is practically any physical question involving the 
explicit knowledge of correlation functions at zero and even more at finite 
temperature. 

The purpose of writing this review was to convince the reader that re- 
search on integrable quantum systems is an attractive and rewarding field 
of science. Hopefully, this goal has been reached and the reader has gained 
from this presentation enough insight to begin or continue his or her own 
investigations. 
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